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Evidence for gravitational quadrupole moment variations in the 
companion of PSR J2051-0827 
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ABSTRACT 

We have conducted radio timing observations of the eclipsing millisecond binary pulsar 
J205 1—0827 with the European Pulsar Timing Array network of telescopes and the Parkes 
radio telescope, spanning over 13 years. The increased data span allows significant measure- 
ments of the orbital eccentricity, e = (6.2 ± 1.3) x 10~ 5 and composite proper motion, 
/i f = 7.3 ± 0.4mas yr _1 . Our timing observations have revealed secular variations of the pro- 
jected semi-major axis of the pulsar orbit which are much more extreme than those previously 
published; and of the orbital period of the system. Investigation of the physical mechanisms 
producing such variations confirm that the variations of the semi-major axis are most prob- 
ably caused by classical spin-orbit coupling in the binary system, while the variations in or- 
bital period are most likely caused by tidal dissipation leading to changes in the gravitational 
quadrupole moment of the companion. 
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1 INTRODUCTION 

PSR J205 1—0827 is the second eclipsing millisecond pulsar system 
discovered after the original black widow pulsar PSR B 1957+ 20 
dFruchter et~al1ll988l) . It was discovered dStappers etaLlll996l) as 
part of a Parkes all-sky survey of the southern sky for low- 
luminosity and millisecond pulsars. The pulsar has a spin period 
of 4.5 ms and inhabits a very compact circular orbit with a very low 
mass companion: m c < 0.1 Mq. The orbital period is Pb = 2.38 h 
and the pulsar and i t s comp anion are separated by just 1.0 Rq. 
Stapp ers et all dl996t 12001 al) investigated the eclipse characteris- 
tics of this system and found that their duration is ~ 10% of the 
orbital period at frequencies below 1 GHz. In addition, assuming a 
mean free electron density of 10 7 cm -3 in the eclipse region, they 
calculated a mass-loss rate of ~ 10~ 14 Mq yr _1 , which is insuffi- 
cient to evaporate the companion completely within a Hubble time. 

Optical observations of t he companion reveale d that it is ir- 
radiated by the pulsar wind dStappersetal.| [l996). Further ob- 
servations and modelling of the variabi lity of the light curve 
(Stap pers et al.ll999l : lsta"ppers et al.l20 01b) determined a binary in- 
clination angle of i ~ 40°, a backside temperature of T ^ 3000 K, 



a companion mass of ~ 0.04 Mq and a radius ~ 0.064 Rq. An 
alternative model with a companion almost filling its Roche lobe 
(~ 0.12 Rq) was also considered and produced almost the same 
results, albeit with a significantly worse fit. For clarification pur- 
poses, a cartoon of the system is given in Figure[T] 

Doroshenko et al. I d200lh presented the most precise timing 
analysis of PSR 12051—0827 using ~ 6yr of radio timing mea- 
surements with the Effelsberg 100-m radio telescope and the 76-m 
Lovell at lodrell Bank. Among the most important measurements is 
the variation of the projected semi-major axis, probably caused by 
Newtonian spin-orbit coupling in this binary system. In addition, 
they measured significant orbital period variations and concluded 
that tho se were created by the same mechanism as in the B 1957 +20 
system ( Applegate & Shaham 1994; Arzoumani an et ai1ll994t) : by 
tidal dissipation of a tidally powered, non-degenerate companion. 

Because black widow (and other eclipsing) systems are sus- 
ceptible to variations of their orbital parameters, monitoring over 
long lengths of time is crucial to understand these sys tems. Tim- 
ing of blac k widow pulsars has be en undertaken by INice et all 
(2000) and lFreireetakl d2003l 120051) , but only very few of these 
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systems have been monitored for a decade or longer. This is of 
particular interest given the large number of black widow pulsars 
recently found i n unidentified Fermi sources IIRansom et a i] |201ll : 
iKeith et al]|201ll . and several discoveries that are soon to be pub- 
lished). 

In this paper we revisit PSR J205 1—0827 with 13 years 
of high-precision timing data and combined datasets from the 
Australia Telescope National Facility's Parkes 64-m radio tele- 
scope and the European Pulsar Timing Array (EPTA) telescopes, 
consisting of the 100-m Effelsberg radio-telescope of the Max- 
Planck-Institute for Radioastronomy, Germany, the 76-m Lovell 
radio-telescope at Jodrell Bank observatory of the University of 
Manchester, UK, the 94-m-equivalent Westerbork Synthesis Ra- 
dio Telescope (WSRT) of ASTRON, the Netherlands, and the 
94-m-equivalent Nancay decimetric Radio Telescope (NRT) of 
the CNRS, France. After describing the properties of our multi- 
telescope data (Section [2}, we present the updated measure- 
ments of the astrometric, spin and binary parameters for the sys- 
tem (Section [3}- Specifically we show the measurement of the 
combined proper motion and of the orbital eccentricity of the 
PSR J205 1—0827 system for the first time and we measure the 
dispersion measure (DM) variations over time. Furthermore, we 
present the extreme variations in the orbital period and projected 
semi-major axis o f the system, which are much larger than the 
ones published by Doroshenko et al. (2001). In Section|4]we rule 
out several possible contributions to the aforementioned variations 
and we discuss the physical mechanisms possibly responsible for 
those. Specifically, we show that gravitational quadrupole coupling 
(GQC) and classical spin-orbit coupling (SOC) can be viable mech- 
anisms for the orbital variations of the PSR J205 1—0827 system, 
under certain assumptions. In light of the new orbital variation 
measurements, we discuss the possible scenarios for the nature of 
the companion and the prospects for observations of the system at 
higher energies, also in Section |4] Finally, in Section [5] we briefly 
summarise our findings. 



2 OBSERVATIONS AND DATA ANALYSIS 

Description of the observing systems and the data acquisition pro- 
cedure of the EPTA telescopes can be found in Lazari dis et al.l 



2009), while for Parkes this information is given in Stappe rs et al.l 



1998). A detailed description on combining pulse ti mes-of-arrival 



(TOA s) from different telesc opes can be found in Ijanssen et al.l 
(2008); Lazaridis et al. (2009); details specific to the data analysis 
used here, are listed below. 

In our timing of PSR J205 1-0827, all TOAs with uncertain- 
ties greater than 20 /is were excluded as they don't impact the 
weighted fit. Furthermore, 573 TOAs that were taken at orbital 
phases between 0.2 and 0.35 were not used because the potential 
excess column density in the eclipsing region may c ause additional 
modulation on the TOAs d erived from these data (Stappe rs et al.l 
ll998l;IStappers et alfcOOlal) . 

Timing was p erformed with the TEMP02 software package 
dHobbs et alj|200fj) . using the DE 405 Solar-System eph emerides 
of the Jet Propulsion Laboratory dStandishl Il998l |2004 and the 
ELL1 binary model dLangeet all 1200 lh . TEMP02 minimises the 
weighted sum of the squared residuals, producing a set of improved 
pulsar parameters and post-fit timing residuals. In order to propa- 
gate unmodelled noise through into the uncertainties of the tim- 
ing model parameters, the TOA uncertainties were multiplied by 
telescope-specific scaling factors to achieve a reduced \ 2 °f unity. 
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Figure 2. Post-fit timing residuals for the last 4.5 years of data. Our timing 
solution is given in TablefJ] The majority of the data are from the Nancay 
radio observatory, but data from Effelsberg, Jodrell Bank and Parkes were 
also used in the preceding nine years, which are not shown. 



These scaling factors all fell between 1.4 and 3.5, where the largest 
factor was applied to the TOAs of the Effelsberg-Berkeley-Pulsar- 
Processor (EBPP). The large reduced x 2 for the EBPP data is most 
likely caused by the fact that this instrument only records total 
power (no polarisation information) and has no RFI mitigation ap- 
plied to it. Arbitrary phase offsets (also called "jumps") were intro- 
duced between data obtained at different frequencies and/or taken 
with different telescopes, to account for frequency-dependent pulse 
shape evolution and observatory-dependent differences in instru- 
mental delays, cable lengths and geodetic position, amongst others. 
These jumps were included as free parameters in the least-squares 
fit performed by TEMP02 and they absorb any non-changing dis- 
persive effects in our data. Therefore we were unable to evaluate the 
averag e DM and hence fixed the value obtained by Stapp ers et al.l 
( 1998) as a constant in our timing model. However, the data does 
contain some sensitivity to time-varying dispersive delays, which 
is fully discussed in Section |33l 

In Table [T] the properties of the data from the different tele- 
scopes are presented. Combination of the TOAs from all the tele- 
scopes provides us with a 13-year dataset with no significant gaps. 
In order to obtain a coherent timing solution for this combined 
dataset, higher order derivatives (up to sixth order) of the orbital 
period and projected semi-major axis must be fitted, leading to high 
correlations between the modelled parameters. These variations in 
orbital period and size are discussed in detail in sections 14.11 and 
14.21 To avoid these correlations and because 85% of our data (and 
all of our highest-precision data) were taken after 2004, we derive 
the timing model, presented in Table [2] based solely on the last 
4.5 years of data, between MJDs 53293 and 54888 (as shown in 
Figure[2}. The remaining older data were exclusively used to inves- 
tigate potential variations in interstellar dispersion (see Section [3~3l 
and to investigate possible causes for the variations in orbital period 
and projected semi-major axis (Section [3.5) , 
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Figure 1. Illustration of the J205 1 -0827 system geometry. The pulsar and companion star are about one solar radius apart, the inclination angle of their orbit 
is approximately 40° and the pulsar wind pushes the companion's wind out, resulting in a bow shock-like denser region of charged particles, which causes 
the system to be eclipsed at given orbital phases and low observational frequencies. Changes in the oblateness of the tidally-locked companion star, possibly 
in combination with alterations in the direction of the spin axis 9, trigger changes in the quadrupole moment that affect the timing of the pulsar. Notice this 
sketch is for illustrative purposes only and is not to scale. 



Table 1. Properties of the individual telescope datasets. 



Properties 


Effelsberg 


Jodrell Bank 


Westerbork 


Nancay 


Parkes 


Number of TOAs 


490 


139 


46 


2954 


51 


Time span (MJD) 


50460-54791 


49989-54853 


54135-54845 


53293-54888 


49982-50343 


Observed frequencies (MHz) 


860, 1400, 2700 


410, 606, 1400 


330, 370, 840, 1380 


1400 


1400, 1700 


Typical integration (min) 


10 


25 


15 


2 


30 



3 TIMING RESULTS 



3.2 Orbital eccentricity 



3.1 Proper motion 

From the data shown in Figure [2] we derived the first significant 
measurement of proper motion in both right ascension and decli- 
nation. Using th e NE2001 model for the Galactic distribution of 
free electrons dCordes & Lazioll2002l) and the pulsar's dispersion 
measure (DM) of 20.745 cm -3 pc, a distance of d ~ 1040 pc has 
been derived. Combining the latter with the composite proper mo- 
tion, fit = 7.3 ± 0.4masyr _1 , the velocity of the system can be 
calculated: 



v t — /-ltd = 36.1 ± 7.5 kms 



(1) 



being consistent with previous estimates ( S tappe rs et al.l 19981) . The 
largest part of the error originates from the DM-derived distance 
where a 20% uncertainty was assumed. 



PSR J205 1—0827 is in a low-eccentricity binary. A precise eccen- 
tricity measurement has been challenging in the past but from our 
timing one has been made for the first time: e = (6.2± 1.3) x 10 . 
This eccentricity value is much larger than expected for such a tight 
binary syste m if plotted in an eccentricity versus orbita l period di- 
agram (as in |Phinnev||l992l : |Phinnev & Kul karni 1994). Although 
the model curves do not extend to very low orbital periods, this 
still constitutes evidence that the PSR J205 1—0827 system at its 
present evolutionary stage experiences processes that do not al- 
low for perfect circularisation. These processes could very well 
be density fluctuations in a convective envelope of the donor star 
(see Section |41~TY Hence the eccentricity of PSR J205 1-0827 
(6.2 x 10 ) is much larger than the expected residual eccentrici- 
ties (^ 10~ 8 ) for binary millis econd pulsars with a similar orbital 
period dLanza & Rodono 2001). 
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Table 2. Timing parameters for PSR J205 1 —0827 for the last epoch of ob- 
servations (MJD range 53293-54888). 



Parameters 


EPTA 


T? 1 writ ntr-PTisinn iT^llflfli 

IVLtLlll *.l >L v. 1 1 > iV 'i 1. LI yj ^\J\.J\J f 


20 h 51 m 07 s 51808(71 


Dprlinatinn X (17i)(M\\ 




11 fmnQvr 




fj,S (masyr -1 ) 


3.2(7) 


v(Hz) 


221.796283737706(1) 


*>(s~ 2 ) 


-6.2639(9) x 10" 16 


P (ms) 


4.50864181828489(2) 


P S -1 } 


1 2733C21X 10~ 20 


IXC1C1C11CC CUUU1 \_1V1J1-'^ 


54091 


MJD range of global timing model 


•J ~J J -J JTQUU 


NJJD range epoch 1 


4QQ eg _ s 1 ^45 


MJD range epoch 2 


51384 - 52332 


MJD range epoch 3 


51967-53117 


MJD range epoch 4 


52790 - 54248 


MJD range epoch 5 


54156 - 54888 


Dispersion measure, DM (cm -3 pc) * 


20.7458(2) 


Orbital period, Pj, (days) 


0.09911024846(2) 


Projected semi-major axis, x (lt-s) 


0.0450725(4) 


7] (= e sin uf) 


5(l)xl0- 5 


K (= e cos ui) 


4(l)xl0~ 5 


Eccentricity, e * 


6(l)xl0~ 5 


Longitude of the periastron, uj (deg) * 


52(12) 


Tasc (MJD) 


54091.0343503(1) 


Pb, epoch 1 (s s" 1 ) 


-1.33(6) x lO" 11 


1 i 1 r\ / .. .. — 1 \ 

-Pb > epoch 2 (s s ) 


0.9(1) X 10 


Pb, epoch 3 (s s _1 ) 


1.8(3) x lO" 11 


Pb, epoch 4 (ss _1 ) 


-1.81(3) x lO" 11 


Pb, epoch 5 (s s" 1 ) 


1.34(6) x 10" 11 


x, epoch 1 (s s" 1 ) 


-1.7(4) x 10~ 13 


x, epoch 2 (s s~ v ) 


-7.9(8) x 10~ 13 


x, epoch 3 (s s~ v ) 


9(1) x 10 -13 


x, epoch 4 (s s~ v ) 


-0.8(2) x 10" 13 


x, epoch 5 (s s" 1 ) 


0.2(2) x 10~ 13 


Solar system ephemeris model 


DE405 


Number of TOAs 


3126 


RMS timing residual (ps) 


12.2 



* The eccentricity and the longitude of the periastron are calculated from 
the Laplace-Lagrange pa rameters, r\ and k. The DM value and its uncer- 
tainty were taken from Stappersetal . (1998) 

Figures in parentheses are the nominal 1 a TEMP02 uncertainties in the 
least-significant digits quoted. 

These parameters were determined with TEMP02, which uses the Interna- 
ti onal Celestial Refe rence System and Barycentric Coordinate Time. Refer 
to lHobbsetai]<2006r) for information on modifying this timing model for 
observing systems that use TEMPO format parameters. 
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Figure 3. DM variations versus time, including the best-fit linear trend. 



phase while keeping the jumps between observing bands and ob- 
servatories fixed. This did not result in significant trends at any 
orbital phase, including egress or ingress-]]. Since this experiment 
required a fit for DM over TOAs that were restricted to a small frac- 
tion of the binary phase, a simultaneous fit to standard orbital pa- 
rameters was impossible. However, the wide frequency range avail- 
able at any given orbital phase implies that even if orbital-phase- 
dependent DM variations did occur, their effect on the TOAs could 
only partially be absorbed by the orbital parameters, implying some 
residual DM variations would still be visible. Simultaneous fitting 
for DM and long-term effects (such as pulsar position and proper 
motion, spin, spin-down and orbital derivatives Pb and x) over lim- 
ited orbital-phase ranges did not affect the DM estimates or the 
other timing parameters, confirming that the long-period terms of 
the timing model are not affected by any potential orbital-phase- 
dependent variations in interstellar dispersion. 

In addition to investigating the DM evolution as a function of 
orbital phase, the DM variations were determined across the entire 
data set. To this end, intervals with adequate multi-frequency data 
were identified by hand and DM values were fitted in each interval, 
while jumps between telescopes and observing bands were kept 
fixed at the values determined from the entire data set. Clearly this 
cannot result in accurate measurements of DM, but it does allow 
variations to be determined precisely. The resulting measurements 
(shown in Fig. [3} demonstrate that no complex DM evolution is 
present in this data set, though a shallow DM decay is measured: 
d(DM) I At = (-4.3 ± 1.4) X 10" 4 cm" 3 yr" 1 pc. This measure- 
ment is barely significant at the 3-cr level. Given the small number 
of high-quality DM estimates that contributed to it and the pos- 
sible underestimation of the DM measurement uncertainties (e.g. 
variations in orbital parameters could have affected the offset be- 
tween non-simultaneous observations), we did not include a time- 
dependent DM variation in our timing model. 



3.3 DM variations 

As described in Section Q] low-frequency observations of this sys- 
tem show eclipses as the pulsar passes through periastron. This may 
imply increased dispersion as a function of orbital phase, which 
could corrupt measurements of binary parameters. To investigate 
this possibility, we took the high-quality, multi-frequency data from 
the final 4.5 years and measured the DM as a function of orbital 



3.4 Companion mass 

Optical observations ( Stap pers et ai1l2001bl) yielded a best fit for 
the orbital inclination angle of about 40°. Depending on the un- 

1 This does not contradict the increased column density in the eclipse re- 
gion but underlines the limited sensitivity of our dataset to DM variations. 
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Figure 4. The masses of the stellar components in the J205 1—0827 system 
as a function of orbital inclination angle. The box indicates the most likely 
parameter space in terms of the neutron star mass and constraints from radio 
eclipses and optical observations (see text). 



known mass of the pulsar this results in a companion mass in the 
range: m c ~ 0.04 — 0.06 Mq (see Figure[4j given the constraints 
from the mass function: 



(m c sin i) 3 
(m p + m c ) 2 



~G 



1.0030 x 10~ 5 M©, (2) 



where the measured observables are the orbital period, Pb and the 
projected semi-major axis of the pulsar orbit, a p sin i. For the rest 
of this paper we shall assume a companion star mass of m c = 
0.05 Mq and consequently a pulsar mass of m p = 1.8 Mq. The 
pulsar mass of PSR J205 1—0827 is likely to be significantly larger 
than the typical neutron star mass (1.35 Mq) obtained from mea- 
surements in double neutron star binaries, given that this system has 
evolved through a lon g (~ Gyr) LMXB phase with sub-Eddington 
mass transfer (see, e.g. Pfahl et al. 2002). Recent work on the mass 
de termination of the original black widow pulsar (PSR B 1957+20) 
bv lvan Kerkwiik et alj J201 lh confirms that the pulsars in these sys- 
tems can accrete significant amounts of matter. For a discussion 
of t he effect of i rradiation on the accretion efficiency in LMXBs, 
see iRitten Q2008). The Roche-lobe radius of the companion star 
in PSR J205 1-0827 changes slightly with the estimated stellar 
masses and is found to be Pl = 0.15 ± 0.01 Rq. However, the 
size (radius) of the irradiated companion star is difficult to deter- 
mine accurately and in the discussion further on we shall assume 
two different values for the filling factor (the ratio of the volume- 
equivalent ^rarHus_of_tlw cornrjanion star to its Roche lobe) of 0.43 
and 0.95 iStap pers et aT]|2001bh . corresponding to stellar radii of 
about 0.064 R Q and 0.14 R respectively. 



3.5 Orbital changes 

In order to monitor their variations over time, values for Pb and 
x were derived for each year of data, where three months of over- 
lap were kept between adjacent years. In doing so, all model pa- 
rameters besides x, Pb and Tasc were held fixed and the timing 
reference epoch was defined to be the centre of each year-long in- 
terval. The fractional changes of these measurements are shown in 
Figure [5] Clearly, five different epochs can be identified, in which 
variations of both Pb and x can be described using only a linear 
trend, as shown in the figure. 
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Figure 5. Changes of the orbital period (top) and projected semi-major axis 
(bottom) versus date. Each point corresponds to close to one year of data, 
with typically three months overlap between adjacent points, though the 
precise amount of overlap varies depending on the density of observations 
and the amplitude of the orbital variations. The horizontal bars indicate the 
intervals over which measurements were made, and the uncertainties of the 
measurements are shown by vertical error bars that are placed at the mean 
TOA of the data contained in the interval. Because of inhomogeneous sam- 
pling, the mean TOA is not necessarily at the middle of the interval. For 
each of the five epochs listed in Table [2] a trend line derived from the re- 
spective Pb an d ^ measurements given in Table f2]is shown. In most cases 
these trend lines fit well through the measurement points, though deviations 
may be caused by rapid variations in orbital parameters, underestimation 
of measurement uncertainties, correlations between parameters and the fact 
that on many of the shorter measurement intervals only part of the tim- 
ing model could be fitted for, incurring potential corruptions in correlated 
parameters. The th ick horizontal bars on the top of the figure indicate the 
timing baseline of Doroshenko et al. (2001), which only really spans the 
first two of these epochs. 



4 DISCUSSION 

4.1 Orbital period variations 

Pb is the observable rate of change of the orbital period and is 
caused by a variety of effects, both intrinsic to the system and 
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caused by kinematic effects relative to the observer. The most im- 
portant contributions are: 



,i GW _i d ■ t 

Pb +Pb +Pb 



+ Pb 



(3) 



As a point of reference, the values for Pb in the two most extreme 
epochs are A = -1.81(3) x lCT 11 and A = 1.8(3) x KT 11 . 

• GW 

The first term, Pb , is the contribution due to gravitational 
wa ve emission. In general relativity, for circular orbits it is given 
by dPetersll 19641) : 



dGW 



1927T 



2n Gm c 
Pb~ 



5/3 



(q+l) 1 / 3 ' 



(4) 



The mass ratio q — 36 has been calculated assuming a pulsar 
mass m p — 1.8 M© and a companion mass m c = 0.05 M© for 
an inclination angle of i = 40°. For PSR J205 1-0827 we find 



Pb 



GW 



-7.5 x 10 14 . This value is about three orders of mag- 



nitude less than the observed value of Pb . 

The second term, Pb , is the Doppler correction, which 
is the combined effect of the proper motion of the system 
(Shklo vsk7i|[T970h and a correction term for the Galactic accel- 
eration. The contribution for the Galactic acceleration at the lo- 
cation of PSR J205 1-0827, P b ° a \ is of order 1.1 x 10 -15 
dLazaridis et alj|2009l) . Using numbers from Table [2] we also cal- 
culate the contribution due to the Shklovskii effect according to the 
following: 



P b Shk = M±MMp b2 , 1 . lxl0 -« 

c 

By summing we yield the Doppler correction: 

P b D = P b Gal + Pb Shk ^2.2xlO- 15 , 



(5) 



(6) 



four orders of magnitude smaller than the measured value. 

An acceleration of the binary system with respect to the Solar 
System Barycentre (SSB) could also be caused by a third massive 
body orbiting the binary system. However, it would affect the or- 
bital period derivative and the spin period derivative in the same 
way. Assuming that the spin period derivative is caused entirely 
by the acceleration, we can estimate the maximal effect this would 



have on P b : (P b /P b ) 



(P/P) 



-3 x 10- 



This is 



seven orders of magnitude smaller than the measured value, thus 
there can be no massive third body orbiting the system. 

The third term, P b , is the contribution from the mass loss 
of the binary. If we consider a circular orbit, no mass loss from 
the pulsa r (m p = 0) and a rate of mass loss from the companion 
m c > (leans 1924), we obtain: 



Pb _ r,m c 
Pb ~ M 
and therefore: 

M P b 

m c = — 

2 P b 



6.2 x 10" 8 M©yr~ 



(7) 



(8) 



where M is the total mass of the system. This numb er is six orders 
of magnitude larger than expected (as calculated by Stapp enTet al.l 
1 1996b and hard to reconcile with the low electron densities mea- 
sured in the eclipse region (Stap pers et al.lll996T) . In addition, for 
the second of the epochs we would need mass injection to the com- 
panion in order to explain the orbital decay - an infeasible scenario 
for this system. 

Tidal forces in tight synchronous binaries may also be respon- 
sible for interactions which would change the orbital period. This 



could be the case if, for example, a companion star suddenly con- 
tracts or expands and thereby changes its spin angular momentum. 
Tidal locking of the orbit would t hen result in exchange of s pin and 
orbital angular momentum, e.g. ITauris & Savoniia (120011) . How- 
ever, in PSR 12051—0827 we do not see reasons for any sudden 
major radial changes in the structure of the companion star - espe- 
cially not changes that would not result in detection of X-ray bursts 
(A. Levine and D. Altamirano, private communication). Hence, we 

• T 

neglect the term Pb . 

Since all the above contributions are much smaller than the ob- 
served variation of the orbital period, we conclude that they must 
originate from the last term of Equation ((3), which represents vari- 
ations of the gravitational quadrupole moment of the companion 
star. 



4.1.1 The Applegate GQC model 

The gravitational quadrupole coupling (G QC) mechanism 
(Applegate 19921: lApplegate & Shahaml 1 19941) has been applied 
successfully to the eclipsi ng binary system PSR B 1957+20 
dApplegate & Sh aham 199 41), expla i ning th e orbital period vari- 
ations seen by Arzoumanian et al. ( 1994). In addition, it was 



proposed by iDoroshenko et al.l ( 120011) as the main mechanism 
for producing the orbital period variations of PSR 12051—0827. 
However, as we will present below, the latter publication con- 
sidered only a small part of these variations, which led to an 

underestimation in their calculations. 

In the GQC model dApplegate & ShaharrJ 1994b . magnetic ac- 
tivity is driven by energy flows in convective layers of the irradi- 
ated companion. This, in combination with wind mass loss, results 
in a torque on its spin, which holds it slightly out of synchronous 
rotation, causing tidal dissipation of energy and heating of the com- 
panion. The resultant time-dependent gravitational quadrupole mo- 
ment (e.g. variations of the oblateness) causes modulation of the 
orbital period on a short, in principle dynamical, time scale. An 
increase in the quadrupole moment causes the two stars to move 
closer together and a decrease in the quadrupole moment results 
in a widening of the orbit - in Applegate's model total orbital an- 
gular momentum is assumed to remain constant and the neutron 
star is treated as a point mass. The details of the hydrodynamic 
dynamo and its activity cycle remain unspecified for our purposes 
(also given the unknown nature of the companion star - see Sec- 
tion |4.4t and it is simply assumed that the variable quadrupole mo- 
ment, AQ, is caused by cyclic spin-up and spin-down of the outer 
layers o f the companion, which lea ds to orbital period changes 
equal to ( Applegat e & Shahamll 1994b : 



AP h 



-9 



AQ 



(9) 



where a is the separation and m c the companion mass. The transfer 
of angular momentum A J to a thin shell of radius R c and mass M 3 , 
rotating with angular velocity Q in the gravitational field of a point 
ma ss m c , will cause the quadru pole moment of the shell to change 



mass m c , will cause the quadrupc 
by dApplegate & Shahamll 19941) : 



AQ 



2 M s Rl 



9 Gm c 



(10) 



where G is the Newtonian gravitational constant and Afi is the 
change in the angular velocity of the shell. In order to produce 
orbital period changes APb you need a variable angular velocity 
given by 
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M s AD. Gm, 
m c S7 



2R? r \R 



,i \ 2 /P h \ 2 AP b 



2tt 



Pb 



(11) 



( Applegate & Sh ahadfl994T) . This procedure produces a luminos- 
ity variation AL given by 



. , tv Gm 2 
AL = — 



3 ^c-fmod V 



Sldr AP b 

SI Pb : 



(12) 



where P mo( j is the period of the orbital period modulation 
and S7dr is the angular velocity of the differential rotation 
( Appleg ate & Shahanj 19941) . 

lApplegate & ShaharrJ dl994l) assume in their model a mass 
M s ~ 0.1m c for the thin shell; an angular velocity Sldr — Afi; 
and luminosity variations at the AL/L ~ 0.1 level. From Fig- 
ure|5]we calculated for PSR J205 1—0827 a total change in the or- 
bital period of APb/Pb — 2.2 x 10~ 7 and a modulation period of 
Pmod — 7.5 yr. As mentioned above th ese two values where u nder- 
and over-estimated, respectively, in Doroshenk o et al.l d200lh be- 
cause of their significantly smaller data set. By assuming an incli- 
nation angle of i = 40 ° and a companion radius of R c ~ 0.064 R© 
( S tapp ers et al.l2001bh for a separation a ~ 7. 1 x 10 8 m and a com- 
panion mass m c = 0.05 Mq, we derive a variation of the angular 
velocity of the companion of: 



AO 

~sT 



4 x 10" 



(13) 



This variation would produce a variable luminosity of AL ~ 1.4 x 
10 33 ergs _1 and according to the model the internal luminosity of 
the companion must be 



L ~ 1.4 X 10 34 ergs~\ 



(14) 



From the optical observations of the companion dStappers et al.l 
2001b) we can estimate a maximum effective temperature of 
T e ff ax ~ 3000 K which gives an internal luminosity of L ~ 
10 30 ergs _1 . We therefore conclude that, under the assumptions 
listed above, the GQC model produces much larger orbital period 
variations than expected and cannot fit the PSR J205 1—0827 sys- 
tem. If the companion was tidally powered, the internal luminosity 
derived from the optical observations would be much higher. One 
reason why GQC does not fit the observations may be the simpli- 
fications of the model itself, i.e. the consideration that all of the 
energy variations appear as luminosity variations without any loss 
(Applegate 1992). Although improvement of the model could pos- 
sibly decrease AL by a factor of a few, the difference we calculate 
is significantly larger. 

The previous calculations can c hange by a large facto r if 
we consider the alternative model that S tap pers et al.l (2001b) pre- 
sented for the companion star, where the latter almost fills its Roche 
lobe (P c = 0.14 R0). Using the values for that model, we derive a 
variation of the angular velocity of the companion of 



AS1 4 
— ~ 7.1 x 10" 4 

and an internal luminosity of 

L~2x 10 31 ergs~\ 



(15) 



(16) 



This number is much closer to the newly calculated internal lu- 
minosity from the optical of L ~ 6 x 10 30 ergs -1 . Thus, under 
specific assumptions the GQC model can explain the orbital period 
variations of PSR J205 1-0827. 



4.2 Changes in the projected semi-major axis 

The observed values of x can be the result of various effects 
dLorimer & Kramej2005l) : 



• GW , ■ D 1 dtA .PM . .m , • Q , • SOC 
Xobs = X + X H — +X + X + X^ + X 

at 



(17) 



and 



For the most extreme epochs, we have: x = —7.9(8) x 10 
x = 9(1) x 10" 13 . 

The first term, i GW arises from orbital shrinkage due to 
gravitational-wave damping. Using Kepler's third law and Equa- 
tion!] 



■ GW 



2 pGW 



-2.6 x 10" 



(18) 



x 3 P b 

dPeterslll964T) . This contribution is much smaller than the current 
measurement. 

The second term, identical to the second term of Equa- 

tion 10. The contribution for the Galactic acceleration is of or- 
der 1.1 x 10~ and the contribution of the Shklovskii effect 
i Shk = x(nl + i4)d/c ~ 6.0 x 10~ 21 . Both of these are very 
small compared to the observed value so this term can be neglected. 

The third term, deA /At, is the contribution of the varying aber- 
ration caused by geodetic precession of the puls ar spin axis and is 
typica lly of order S7 gcod P/P b « 7.4 x 10" 17 dDamour & TavloJ 
I1992I) . For a recycled pulsar, like PSR J205 1-0827, the spin is 
expected to be close to parallel to the orbital angular momentum, 
which further suppresses this effect. Hence, the contribution is at 
least three orders of magnitude smaller than the observed value. 

The fourth term, i PM , represents a variation of x caused by a 
change of the orbital in clination while the binary system is movin g 
relatively to the SSB dArzoumanian et alj[l99tt lKopeikinlll996h . 
This effect is quantified by the following equation: 



i PM = 1.54 x 10" 



1 x cot i(—n a sin fi asc + fj,$ cos Sl a 



(19) 



where S7 asc is the position angle of the ascending node. The quanti- 
ties x, )j, a and jj,g are expressed in seconds and milliarcseconds per 
year, respectively. The maximal contribution of the proper motion 



i™ = 1.54 x 10" 



' x (ji 2 a + jis) 1 coti. 



(20) 



For an inclination angle of % = 40°, we get £™ x ~ 6.0 x 10~ 17 . 
Thus, this term is also very small compared to the measured x. 

The fifth term, x m , represents a change in the size of the orbit 
caused by mass loss from the binary system. Using Equation (|7} 
and Kepler's third law, we calculate the rate of mass loss from the 
companion to be seven to eight orders of magnitude larger than 
expected. 

The sixth term, if 1 , is caused directly by the change in 
quadrupole moment as described in Section |4~T1 Because orbital 
angular momentum is conserved in the Applegate model, a change 
in the orbital period is related to a change in the size of the orbit by 



Aa 



.AP b 



(21) 



It can be seen from Figure [5] that this contribution falls short by 
four orders of magnitude to explain the observed changes in the 
projected semi-major axis x. 

As all other contributions are much smaller than the observed 
variation of the projected semi-major axis, they must originate from 
the last term of Equation dl7t : the classical spin-orbit coupling 
(SOC) term. 
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Table 3. Quadrupole moment variations for the five epochs, as derived from 
the Pt, variations, alongside the x variations as observed. 



le-06 



Epoch 


AQ 

(f0 45 gcm 2 ) 


n h AQ 
(lCT 11 s" 1 ) 


x/x 

(lcr 11 s^ 1 ) 


1 


11.2 


7.7 


-0.38 


2 


-4.6 


-3.2 


-1.8 


3 


-11.3 


-7.7 


2.0 


4 


14.3 


9.8 


-0.18 


5 


-5.3 


-3.6 


0.044 



4.2.1 The Applegate GQC model with SOC 

Through SOC, the quadrupole of a rapidly rotating companion 
leads to apsidal motion and precession of the binary orbit. This in 
turn causes a variation of the longitude of periastron (which is im- 
possible to measure in a system with such small eccentricity) and 
of the projected semi-major axis, according to: 

x SOC = xnbQ cot i sin 9 cos 9 sin $ , (22) 

( Smarr & Blandford 1976; Wex 1998), where rib = 2n / Pb is the 
orbital frequency, 9 is the angle between the spin and orbital angu- 
lar momentum and $ is the longitude of the ascending node with 
respect to the invariable plane (plane perpendicular to the total an- 
gular momentum). The dimensionless quadrupole Q is related to 
the quadrupole Q by 



2 V a 



•h = 



3Q 
m c Rr 



(23) 



with J2 the dimensionless measure of the quadrupole moment. 
Since no uniform change in x is detected, we can consider the 
companion to consist of a long-term stable component with spin 
axis aligned to the orbital angular momentum axis (i.e. 9 ~ 0) and 
in addition to that, a part of the star (like e.g. an outer shell) that 
changes its quadrupole moment by AQ and has an effective and 
variable angle 9 independent of the rest of the star. This part will 
then give rise to a variation x/x of the order ribAQ. 

In Section |4.1.1| we have described how changes in the gravi- 
tational quadrupole moment of the companion could cause the ob- 
served Pb variations. Even though these direct GQC effects are in- 
sufficient to explain the x variations (i.e. = as explained 
above), through the SOC mechanism they might have a more sub- 
stantial effect. 

Extending Equation|9]to arbitrary 9 and using Equation|23] we 
have: 



AP b 



-2AQ 1 



3 . 2 Q 
— sin 9 

2 



(24) 



Combining this with the observed Pb values (Table |2j and assum- 
ing 9 = for now, we get the maximum AQ values presented in 
Table [3] Comparison to the x/x values from Figure [5] also listed 
in Table [3] shows that while the order of magnitude is rougly cor- 
rect, the sign does neither correlate nor anti-correlate, which means 
that the orientation of the quadrupole-moment changes (i.e. 9) must 
vary strongly from epoch to epoch, which may be unphysical. 

In order to examine more closely if the orbital variations can 
be produced by GQC and SOC arising from the same AQ, we have 
plotted in Figure [6] how |AQ| depends on 9, assuming $ = 45° 
and using i = 40°. The curves are plotted for SOC, as derived 
from Equation[22] and for GQC, as derived from Equationf24] The 
intersection of the two lines clearly shows that for several values of 



^ le-07 



le-08 



le-06 




45 90 135 180 
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Figure 6. Quadrupole moment changes versus 6 as produced by SOC (solid 
line) and GQC (dashed line) for epochs two (top) and three (bottom). 



9 a AQ could explain both the orbital period and projected semi- 
major axis variations for PSR J205 1 —0827, though these 9 values 
change significantly with epoch and we do not see how this easily 
fits with the phy sical mechanism behind the thin shell model of 
Applegatel \ 19921) (see also our Section E.l.lt . 



4.2.2 Spin precession of the companion star 



In Section 14.2.1 1 we assumed a mostly stable companion star with 
an outer shell responsible for the changing quadrupole moment. 
An alternative explanation for the variations in Pb and x could be 
provided by limited quadrupole changes (as needed for the Pb va- 
rations, see Section |4. Lit combined with an overall small tilt 9 of 
the entire companion star with respect to the orbit. This tilt causes 
a precession of both the star and the orbit and it allows the entire 
quadrupole moment of the star Q to cause variations in x according 
to Equation [22] We stress that in this scenario there need not be a 
physical mechanism relating the spin precession (which causes the 
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Table 4. Quadrupole changes (orientation and magnitude) that explain the 
changes in Pb and x for the five epochs. 



Epoch 


0(deg) 


AQ/Q 


1 


-0.14 


+0.041 


2 


-0.64 


-0.017 


3 


+0.73 


-0.041 


4 


-0.06 


+0.052 


5 


+0.02 


-0.019 



x variations) to the quadrupole moment changes (that induce the 
Pb changes^- 

To evaluate this model, we estimate the quadrupole moment 
of the co mpanion star assuming a non-degenerate structure^ and 
following lLanza & Rodonol fl999l) : 



5 R? [T-(l-3y)M] 
V 6£ GM 



(25) 



where R and M are the radius and mass of the star, T is its rota- 
tional kinetic energy (T = IQ 2 /2), M is its total magnetic energy, 
r\ is the fraction of this magnetic energy in a field directed along the 
rotation axis and £ is a dimensionless parameter depending on the 
density stratification and the spatial distribution of the perturbing 
forces in the convective layer of the star (it can be calculated from 
the gyration radius and the apsidal motion constant using stellar 
models). Here we shall use £ ~ 5. The rotational frequency of the 
outer layers of the star is taken to be that of the orbit (caused by 
synchronisation) and in Applegate's model M = 0. Adopting a 
companion star radius R — 0.064 R© yields: 

Q = 3.4 x 10 47 gcm 2 . 

Consequently, the quadrupole moment Q of the star is one to two 
orders of magnitude larger than the variations required for the ob- 
served Pb changes. Using this value in Equation [22] results in the 
angles 9 required for the observed x variations. These values can 
in turn be used to correct the AQ values shown in Table[3]for the 
angular dependence, though this effect is all but negligible. The re- 
sulting changes required in the tilt and quadrupole moment of the 
companion star are collated in Table [4] 

We conclude that with small changes in the orientation and 
magnitude of the quadrupole moment, we can explain all observed 
parameter changes. The change in the sign of x from one epoch to 
the other can be easily explained by a small oscillation of the sym- 
metry axis of the quadrupole. Moreover, if we use an alternative 
scenario with a larger filling factor for the radius of the compan- 
ion star, these numbers will get even smaller since Q will be larger. 
We do not, however, know of any physical mechanism that would 
cause the required tilt of the star or its rapid changes. 



2 Note the changes in 9 will have some contribution to the Pb changes as 
well. Based on Equation 5 of Appl egatel j 19921) we determined this effect 
and for all epochs save the third, this effect is at or well below the few 
percent level of the Pb changes induced by A Q. In the third epoch, which 
requires the largest angle 8, it contributes 1% of the observed Pb variations. 

3 Deviations of a semi-degenerate nature should not be important in this 
case. 



4.3 Alternative Models for the Gravitational 
Quadrupole-Moment Variations 

The Applegate model is one of a number that have been applied to 
explain the cyclic modulation in the orbital period of magnetically 
active close binaries. In Applegate's model a rather large fraction 
of the stellar luminosity is required for its operation. Hence, one 
might be able to detect changes in the stellar luminosity in phase 
with the orbital period modulation. The model of Lanza & Rodono 
(e.g. lLanza et aL 1 11998k lLanza & Rodonol Il999h considers varia- 
tions in the azimuthal B-field to explain the variations in oblate- 
ness. A change in the azimuthal field intensity can also produce a 
change in the quadrupole moment by changing the effective cen- 
trifugal acceleration. However, one should keep in mind that it re- 
mains to be established what kind of companion star we have in 
PSR J205 1-0827. 



4.4 The Nature of the Companion Star 

There seem to be three possibilities for the nature of the companion 
star: 

I a white dwarf (WD), 
II a brown-dwarf-like star and 
III a semi-degenerate helium star (He). 

In the following we briefly discuss each of these possibilities. 

For a WD companion star of M c ^ 0.04 Mq, the mass- 
radius relation of a non-r elativistic degenerate Fermi-gas (e.g. 
Shapi ro & Teukolsky|[T98l yields an upper limit for the radius of 
R c ^ 0.025 Rq. This value is much less than what is estimated 
in optical since a combination of the minimum value for the filling 
factor of 0.43 and the Roche-lobe size of Pl = 0.15 Rq, yields 
a minimum optical radius of P c ^ 0.064 Rq, However, because 
of irradiation by the pulsar wind and tidal dissipation of energy in 
the WD envelope it is likely that the WD is bloated in size. The 
effect of a thermally bloated WD in a close binary system has re- 
cently been detect ed in a transiting source by the Kepler satellite 
dCarter et al.ll201 lb . In this case the young (few hundred Myr) hot 
white dwarf is bloated by a factor of seven in size. The companion 
star in PSR J205 1 —0827 is believed to be much older, however it is 
quite possible that an extended H-rich atmosphere may exist in this 
star from the effects mentioned above. On the other hand, the actual 
effect of the pulsar wind remains uncertain and may not be as effi- 
cient a s expected previously. For example, Equation (2) of iTavanj 
dl992l) for the irradiation of the companion as induced by the pul- 
sar wind, yields a surface temperature of the compan ion of 6200 K. 
Howev er, the observed temperature is ^ 3000 K (Stapp ers et alj 
l2001bl) . 

Brown dwarf mod els have been applied 
(Bildsten & Cha krabartvl |200lh to explain the nature of the 
~ 0.05 Mq donor star in the accreting millisecond pulsar system 
SAX J1808.4-3658. Brown dwarfs have a relation between mass 
and radius that is acceptable from observational constraints, also 
for PSR J205 1—0827. However, the obtain ed mass-radius r elatio n 
was based on brown-dwarf models by Chabrier et al.l fcOOOl) . 
calculated for isolated low-mass stars. We find it questionable 
if the present core remnant in PSR J205 1—0827 can simply be 
described by the equivalent of a brown dwarf. The original ZAMS 
progenitor system of PSR J205 1—0827 definitely did not consist 
of a B-star and a brown dwarf: the initial sub-stellar mass of a 
brown dwarf is < 0.08 Mq (the limit for H-ignition) and the 
progenitor of the neutron star had a mass > 10 Mq. Hence, the 
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initial binary would have had a mass ratio of < 1/100, which is 
not only disfavoured by formation, but also unlikely given that 
this binary should undergo a common envelope phase followed 
by survival of a supernova explosion and finally spin up the 
pulsar with so little mass. A more reasonable possibility is that 
PSR J205 1—0827 (and many other eclipsing binary millisecond 
pulsars) is the outcome of a short orbital period, low-mass X-ray 
binary with un-evolved main se quence stars. Such s y stem s have 
been studied in detail by, e.g., iPvlvser & Savoniiel dl988l) and 



Ergm a~et"al1 dl997t) . The latter authors also discussed the final fate 
of a system like PSR J205 1-0827: an ultra-compact LMXB where 
the companion fills its Roche lobe (again), because of orbital 
angular momentum losses, and eventually may undergo tidal 
disruption leaving behind a single or planetary millisecond pulsar 
(e.g. a system like PSR B1257+12. IWolszczan & Frailll992h . 

Finally, it should be noted that the companion star in 
PSR J205 1—0827 could be closely related to the semi-degenerate 
helium star companions in ultra-compact A M CVn systems (e.g. 
Ilben & Tutukovlll99ll ; lNelemans et alj|200lh. One ex ample is the 
system SDSS J0926+3624 dCopperwheat et al.ll2010l) which has a 
similar companion mass of 0.035 M© (Pb = 28min). Although 
the accreting compact object in AM CVn systems is a WD, a simi- 
lar evolution is expec ted for systems with an accreting neutron star 
dSavoniie et aljfl986l) . Also, in this case the expected mass-radius 
relations yield acceptable values for the radius in PSR J205 1 —0827 
- especially if one accepts the possibility that a pulsar wind and 
tidal dissipation could lead to a somewhat bloated radius. 



4.5 High-Energy Emission Prospects for PSR J2051 0827 

The detection of X-rays from PSR J205 1—0827 could yield vital 
information about its orbital evolution, the state of the compan- 
ion star and interactions between the pulsar wind and the evapo- 
rated material of the companion star. However, no X-rays are de- 
tected from PSR J205 1 -0827 in the RXTE All-Sky-Survey (A.M. 
Levine, private communication). The RXTE ASM sensitivity in 
the 2-10 keV band is ~ 30mCrab, corresponding to an energy 
flux of F x ~ 7 x 10~ 10 erg s -1 cm -2 . Given the distance to 
PSR J2051— 0827 (d ~ 1.04 kpc), we can estimate an upper limit 
to the mass accretion rate: 



4tt gT F x R d 



(26) 



Assuming a neutron star mass m p = 1.8 Mq and radius R p = 
10 6 cm, we obtain a conservative upper limit for the potential ac- 
cretion rate m < 6.0 x 10~ 12 Mq yr _1 . Keeping in mind that the 
expected mass transfer rate for a Roche-lobe-filling ~ 0.05 Mq 
companion star is expected to be much lower than this limit - and 
given the fact that we do observe a radio pulsar - it is not a sur- 
prise that PSR J205 1 -0827 is not detected in the RXTE ASM data. 
Chandra and/or XMM data would impose far more stringent con- 
straints; or might even provide a detection. 

One must bear in mind that the lack of detected X-rays does 
not exclude mass loss from the companion star. It is possible that 
the pulsar wind (which is probably enhanced in the equatorial re- 
gion towards the companion star) is able to prevent any accretion 
onto the neutron star. However, in this scenario it needs to be in- 
vestigated thoroughly whether or not a shock front would lead to 
acceler ation of protons and subsequent production of high-energy 
7-rays ( Hardi ng &Gaisserlll990t IStappers et afll2003h . 

Considering 7-ray emission from this system, with a distance 
of ~ 1.04kpc and a spin-down luminosity E = 4tt 2 IP/P 3 = 



5.5 x 10 33 erg/s, this pulsar is a good candidate for detection wit h 
the Large Area Telescope of the Fermi satellite dAbdo et alj|2009h . 
However none of the Fermi sources in the fi rst-year catalogue lie 
within a radius of 3° of PSR J2051-0827 dAbdo et alj|20ic|) . It 
would therefore be interesting to see if PSR J205 1—0827 is de- 
tected by Fermi with accumulated data. 



5 SUMMARY 

We have presented a timing update on PSR J205 1—0827 and shown 
that the variations of the orbital period and projected semi-major 
axis of the binary system are far more extreme than described in 
earlier work. We have analysed all possible causes for these vari- 
ations and found that a combination of gravitational quadrupole 
coupling and spin-orbit coupling are the most likely origins of this 
behaviour, though accurate modelling efforts are still required. We 
furthermore discussed the nature of the companion star and con- 
clude that a semi-degenerate helium star is more likely than the 
white-dwarf alternative and, finally, we expect high-energy 7-rays 
may be visible in this system if a pulsar wind prevents accretion 
onto the neutron star - which in turn is expected based on the lack 
of X-rays. 
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